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. Abstract 

OC • 

CN ' We present analytical results on the two-loop anomalous dimensions of currents for baryons, 

containing two heavy quarks J = [Q 17 CTtQ^]T ' q k '£ijk with arbitrary Dirac matrices T and T' in 
the framework of NRQCD in the leading order over both the relative velocity of heavy quarks 
and the inverse heavy quark mass. It is shown, that in this approximation the anomalous 
■ dimensions do not depend on the Dirac structure of the current under consideration. 
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1 Introduction 



The necessary feature of QCD applications to various fields of particle physics is a study of a scale 
dependence for operators as it is governed by the renormalizat ion-group (RG). In the present paper 
we investigate the RG properties of currents for baryons with two heavy quarks in the framework 
of Non-Relativistic Quantum Chromodynamics (NRQCD) ||1|,[01 an d its dimensionally regularized 
version ||. In the two- loop approximation we analytically calculate the anomalous dimensions of 
currents associated with the ground-state baryons, containing two heavy quarksQ. The dependence 
of QCD operators and matrix elements on the relative velocity v of heavy quarks inside the hadron 
as well as on the inverse heavy quark mass 1/Mq can be systematically treated in the framework of 
justified effective expansions in QCD. So, we apply the expansion in 1/Mq, as it was developed in 
Heavy Quark Effective Theory (HQET) J7], |8], EJ for operators, corresponding to the interaction of 
heavy quarks with the light quark. For the heavy-heavy subsystem, the power tool is the NRQCD- 
expansion in both the relative velocity and the inverse mass. Here we consider the leading order in 
both v and 1/Mq, that can serve as a good approximation for the anomalous dimensions of currents 
under consideration. 

The anomalous dimensions of composite operators can be desirably used in QCD sum rules [TU|. 



which will allow us to evaluate the masses of these baryons together with their residues in terms 
of basic non-perturbative QCD parameters. For example, calculating the two-point correlators of 
baryonic currents in the Operator Product Expansion (OPE) in NRQCD, we have to insert the 
anomalous dimensions, obtained here in the static approximation, to relate the result to QCD. 

This procedure is caused by a different ultraviolate behaviuor of loop corrections in the full QCD 
and the effective theory. The latter contains the divergences absent in QCD, since it was constructed 
in the way to provide correct infrared properties of local QCD fields. The regularized quantities of 
effective theory depend on the normalization point under the RG equations with the corresponding 
anomalous dimensions. The ambiguity in the initial conditions of such differential equations is 
eliminated by the matching to full QCD at a scale, which is generally chosen as the heavy quark 
mass. The latter procedure means, that using the effective theory, we can systematically take into 
account the virtualities greater than the heavy quark mass. So, the knowledge of the two-loop 
anomalous dimensions is also important, when one discusses the matching of baryonic currents, 
obtained in this approximation with the corresponding currents in full QCD. 

Our analysis in this paper is close to what was presented in [|ll|], devoted to the baryons with a 
single heavy quarkP]. While being very similar, these analyses also have some differences, which we 
would like to stress. The main technical obstacle of calculations is related to that the kinetic term 
is thought to be a necessary ingredient in the quark propagator for the evaluation of RG quantities 
in NRQCD, unlike to HQET, 
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If a hard cut-off is used (/x <C m), we can easily see that such NRQCD-calculations can be per- 
formed just like in HQET, since k° 3> k 2 /m in the ultraviolet regime. However, if the dimensional 
regularization is used, the high energy modes (k > m) are not explicitly suppressed and they give 
non-vanishing contributions. This can be seen because the behavior of the NRQCD propagator 

1 We do not consider the problems concerning the spectroscopy, decays and production mechanisms of baryons 
with two heavy quarks. This can be found in 0,S and correspondingly. 
2 We generally accept a set of basic notations used in [IE. 
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changes at energies greater than the mass. In spite of this, one would like to use dimensional regu- 
larization because it keeps all of the QCD symmetries and, moreover, the calculations are technically 
simpler. 

The difference between NRQCD and HQET can be explicitely highlighted in the consideration 
of effective Lagrangian, derived to the tree level in the 1/m-expansion: 

£nrqcd = ^(iD + £-)il> + -L^D^-^^tr-B^ 
\ 2m J 8m A 2m 

- -^r ift (D ■ E — E ■ D)ib — -^r V>V ■ (D x E - E x D) ib 

8m 2 8m 2 

+0(l/m 3 ) + antiquark terms + £ii g ht (2) 

For a single heavy quark, interacting at low virtualities D ~ ^qcd, the kinetic term is suppressed 
and can be treated perturbatively. This results in the HQET prescription to the heavy quark 
propagator. However, in the heavy-heavy system there is the Coulomb-like interaction, wherein 
D° ~ D 2 /m ~ a 2 s m. Therefore, we must include the kinetic term into the initial "free" Lagrangian 
of NRQCD. So, the loop corrections in a s look to be different in HQET and NRQCD. Nevertheless, 
the physical reason to distinguish these effective theories is still the Coulomb-like corrections near 
the production threshold, which should make no influence on the ultraviolate properties. We would 
note that the question is, in a sense, analogous to that in the theory of massive gauge fields in the 
spontaneously broken theories, where the explicite introduction of mass seems to destroy good RG 
properties of massless vector fields (the question was removed by the appropriate redefinitions of 
fields due to the surviving the gauge invariance). 

Several authors have addressed the similar problem of NRQCD in the connection to matching 
calculations WM, and recently an appealing solution has been proposed [Of: it is claimed that the 



matching in NRQCD using the dimensional regularization should be performed just like in HQET, 
namely the kinetic term must be treated as a perturbation vertex: 

Ilk 2 

T- + 7r-7T^ + - (3) 



The derivation is based on the appropriate redefinition of the heavy quark field fl3| : 



Q _^ {1 Dj g^G^ | Djjjv ■ D)Dg± | gvxD±*G aX 
8m 2 16m 2 16m 3 16m 3 
. cr al3 Dl{iv D)Dj _ . gv x <7 aP DlG^ 

16m 3 16m 3 ' 

where the a matrices are projected by P v aP v , P v = and = — v^v ■ D. The substitution 
converts the HQET Lagrangian to the NRQCD one, so that the loop renormalization of perturbative 
terms is the same. 

Here, we propose to use the same prescription for the heavy quark propagator as it stands in 
(|3]), not only in the matching procedure, but also for the calculations of anomalous dimensions 
for the NRQCD currents in MS-renormalization scheme. To support this point let us consider the 
matching procedure in some details. The matching condition can be written down as 

y — 1 yon— shell yh.m. tV yon— shell y — 1 p(0) tr\ 

^ 'J,QCD^2,QOD Zj V,QCD l QCD — { ~'0^2,NRQCD Zj J,NRQCD l NRQCDi \°) 
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yh.m p' p(0) in\ 

^V&CD 1 QCD — 1 QCDi V°) 

where Zy m ' denotes poles, associated with the hard momenta region for the bare single-particle 
irreducible vertex Fq CD in full QCD, Zjq CD and Z JNR q C d are the renormalization constants of 
currents in QCD and NRQCD, correspondingly, Z 2) qcd and Z 2 ,nrqcd include the renormalization 
of wave functions, and, finally, ^^rqcd denotes the bare vertex in NRQCD. On this stage we use 
prescription (|3|) for the treating the heavy quark propagators. On the other hand, one can write 
the following indentity 

p y — 1 7MS yh.m. ys.ra. p" (rj\ 

QCD — A J,QCD Zj 2,QCD Zj V,QCD Zj V,QCD l QCDi \< J 

where we have collected all divergences in Z-factors and use the convention of (|3D for the expansion 
of heavy quark propagators in powers of the kinetic term. Zy 1 Q , CD denotes the contribution from a 
small momenta region. Calculating the contribution from the small momenta, we have to set the 
external legs to be off-shell, in order to exclude the contribution from the infrared region as it was 
done in the case of matching. To proceed further, let us introduce the following definitions 

yon— shell 7MS y ( o\ 

^2,QCD ~ Zj 2,QCD Zj inf.r., \°) 

yon— shell yMS y (q\ 

Zj 2,NRQCD ~ Zj 2,NRQCD Zj inf.r. } \V) 

where Z in f, r , is a contribution to the wave-function renormalization from the infrared region, which 
is the same in both theories. Using these notations and the fact, that Z^nrqcd = 1 ; we can rewrite 
Eq. (0) as 

p y — 1 yon— shell yh.m. ys.m. 7MS p" /1 r\\ 

L QCD — A j,QCD Zj 2,QCD Zj V,QCD Zj V,QCD Zj 2,NRQCD l QCD- \ LU ) 

Now we can easily see from Eqs. ([5]) and (TU), that the NRQCD anomalous dimensions in the 



MS-renormalization scheme can be computed in two ways: either from the matching condition 
(|5|) or using the HQET Feynman rules and setting the external legs off-shell in order to avoid 
the infrared divergencies. We have explicitly checked this conjecture to one-loop for the heavy- 
heavy vector current, however, for a full confidence we feel a need for such a check in the two-loop 
approximation. 

So, in our approach we exploit the same reasoning for the calculation of RG quantities and 
work in the leading order of this expansion. Moreover, it is theoretically sound, because in the MS- 
renormalization scheme used, the anomalous dimensions of currents do not depend on the masses 
of particles. The following fact also supports our claim: the values of Wilson coefficients, calculated 
in the matching procedure, are directly connected to the anomalous dimensions of operators mul- 
tiplying these coefficients in the Lagrangian. And, finally, the high energy behavior in the effective 
theory with several scales does not depend on relative weight of the lower scales. Thus, we only 
need 

m> \p\,E,A QC d, (11) 

where there is no matter what are relations between |p|, E and Aqcd- 

So, in our calculations we use the HQET propagators for the heavy quarks, setting the quark 
momenta in a way to avoid infrared divergencies. As will be explained in details below, the two-loop 
contribution to the anomalous dimensions of currents under consideration consists of three parts. 
The first corresponds to the set of graphs, wherein the two-loop contributions are associated with 
one of the heavy- light subsystems. For this contribution we use the result of [jTTJ . Then there is 



the subset of two-loop graphs that are associated with the heavy-heavy system. The expression 
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for this contribution is a generalization of what was obtained in |14]]. And, finally, there are the 
irreducible contributions, where the two-loops connect all three quark lines. This contribution is 
calculated in this paper. We evaluate the two loop diagrams with the use of package, written by us 
on MATHEMATICA, and the recurrence-relations in HQET [15]. 

This paper is organized as follows. In section 2 we discuss the choice of currents for the baryons 
with two heavy quarks and give some comments on the renormalization properties of composite 
operators under consideration. In section 3 we furnish some remarks on the anomalous dimensions 
and present the results on the one-loop anomalous dimensions. In section 4 we discuss general 
features of two-loop renormalization procedure and present our two-loop results. We work in the 
MS-renormalization scheme throughout the paper. As concerns the treatment of 75 we will show 
that the final expression does not depend on the scheme used. Section 5 contains our conclusion. 



2 Baryonic Currents 

The currents of baryons with two heavy quarks S* c , S£ 6 and Hj°, where means different charges 
depending on the light quark charge, are associated with the spin-parity quantum numbers j p = 1 + 
and j p = + for the heavy diquark system with the symmetric and antisymmetric flavor structure, 
respectively. Adding the light quark to the heavy quark system, one obtains j p = \ + for the S' b * 
baryons and the pair of degenerate states j p = | + and j p = | + for the baryons S* c , E bc , E£ 6 and 
, 'EH, — s>& • The structure of baryon currents with two heavy quarks is generally chosen as 

J = i)' 1 ( T? () J V <i' ,js.-- (12) 

Here T means transposition, C is the charge conjugation matrix with the properties C^ T C~ l = — 7^ 
and C^JC^ 1 = 75, i,j, k are colour indices and r is a matrix in the flavor space. The effective static 
field of the heavy quark is denoted by Q. To obtain the corresponding NRQCD currents one has to 
perform the above-mentioned redefinition of local field. But as we are working in the leading order 
over both the relative velocity of heavy quarks and their inverse masses, this local redefinition does 
not change the structure of the currents. 

Here, unlike the case of baryons with a single heavy quark, there is the only independent current 
component J for each of the ground state baryon currents. They equal 



J^ Q , = [Q iT CT l5 Qi'}q k e l3kj 

Jz« QQ = [Q lT Cr^]- ll5 q k e ljkl (13) 
Js* Q% = [Q iT CT 1 Qq q % k + p[Q iT C 7 Qi]-~fq k e ijk , 
where Je*° q satisfies the spin-3/2 condition "yJs^ Q = 0. The flavor matrix r is antisymmetric for 5 



be 

and symmetric for Sqq and The currents written down in Eq. (6) are taken in the rest frame 
of hadrons. The corresponding expressions in a general frame moving with a velocity four-vector 
v M can be obtained by the substitution of 7 — > 7^ = 7^ — vv^. 

Now we would like to give some comments on the renormalization properties of these currents. 
As we have the only one light leg in this problem, all of 7 matrices, which will appear in calculations, 
will stay on a single side of our composite operators, not touching their Dirac structure. This will 
lead to the fact that the anomalous dimensions of all our currents in this approximation are the 
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same, i.e. they do not depend on T-matrices in (4). From this reasoning, we also can conclude that 
the result does not depend on the 75 scheme used. 



3 Common notations in renormalization 

The local operators Oo composed of bare physical fields contain the ultra-violet divergences, which 
can be absorbed by the renormalization factors Zq, being a series in powers of coupling constant, 
so that O = ZqOq is a finite quantity, while the regularization parameters do not tend to peculiar 
values. In the dimensional regularization using the MS-scheme of subtractions in D = 4 — 2e 
dimensions Zq is expanded in inverse powers of e, so that 

oo m 00 1 

m=l k=l k=l 

The dependence on the dimensionful subtraction point /x defines the anomalous dimension of renor- 
malized operator O 

dla.Z(a(n),a;e) 

7= — dW) — ' (15) 

where a is the renormalized gauge parameter in the general covariant gauge (with a gluon propagator 
proportional to —g^ u + (1 — ajk^ky/k 2 ) and a(/x) is the renormalized coupling constant in four- 
dimensional space, so that 

a = a(fj,)fi 2e Z a (a(fi), a; e), a = aZ 3 (a(/x), a; e), (16) 

and the corresponding Z^^-factors determine the anomalous dimensions, which are generally de- 
noted by {— P, —5}, respectively. 

The 7-quantities are finite at D — ► 4, so we define the coefficients of series 



00 



7<m ' (17) 

m=l 

One can check that [17| 

7 = "2^-, (18) 

and for k > 

_ 2 dZ k+1 = (_ g 9 _ § 9 \ Zk / 1Q N 
d\na s \ d\na s d\na) 

The latter provides the consistency condition, when the former produces a simple extraction of the 
anomalous dimensions to the two-loop accuracy 

7 (1) = -2Z M and 7 {2) = -4Z 2]1 . (20) 
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3.1 One-loop result 



Consider the one-loop renormalization of currents of baryons with two heavy quarks. In the MS- 
scheme with D = 4 — 2e space-dimensions we have the following squares of renormalization factors 
for the bare quark fields: 

Z « =1 -°°!!? Z ^ 1 + < - 3 -^WFe' (21) 

where we use the usual definitions for SU(N), i.e. C F = (iV c 2 -l)/2iV c , C A = N c , C B = (N C +1)/2N C , 
and Tp = 1/2 for N c = 3, Np being the number of light quarks. One- loop MS-results for the factors 



Z a and Z 3 have been given e.g. in [17 



1 + 
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—C A - —TpNp 
3 3 
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^C A - -TpNp 
6 3 



(22) 
(23) 



The bare current is renormalized by the factor Zf. 

J = (Q^CTrQ )r qo = ZgZ^ZyJ = ZjJ, (24) 
which straightforwardly means that 

7 J = 27q + l q + Jv, (25) 

i.e. the anomalous dimension of the full current J is a sum of three terms given by the renormal- 
ization of the light and heavy quark fields, and the renormalization of the vertex. 
For the vertex {Q1CYtQq)Y' q Ql we find 



which results in 



'V 



7? } 



1 + 



4vre 



■(3a — 3), 



-2C B (3a-3). 

The one-loop anomalous dimensions 7^ and 7q are equal to 

T « = C F a, 7 « = C F (a - 3). 
Thus, the one-loop anomalous dimension of the baryonic current is given by 



«.s- 



47T V 



( - 2C B (3a - 3) + 3C F (a - 2)) + 0(a 2 s ). 



(26) 
(27) 

(28) 

(29) 



4 Two-loop calculations 

In this section we apply the two-loop renormalization of the baryon current with two heavy quarks in 
the MS-scheme and restrict ourselves by the Feynman gauge. The two-loop anomalous dimensions 
of the quark fields are given by H, g]|, H || . 

7? = C F {^Ca - 2T F Np - \ C ^j , = C F (-^C A + f T pNp^ . (30) 

Since the baryonic currents are renormalized multiplicatively in the effective theory^, the Dirac 
3 See discussion in ref. |TT| . 
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structure of vertex repeats the Born-term. Technically we perform the calculations in terms of 
bare coupling and gauge parameter, so that to isolate the two-loop contribution to the anomalous 
dimension, we need also the one- loop result, wherein we have to include the one- loop expressions 
written down through the renormalized quantities a s and a, which will add the contribution to the 
corresponding a^/e-term. The procedure described leads to the relations: 

Z\,l — ^1,1) ^2,2 = ^2,2) -^2,1 = ^2,1 — V%, 1^1,0- (31) 

As expected, Z 2) i has to include the one-loop contributions. 

In the Introduction we have described three subgroups of two-loop diagrams for the vertex, 
which can be expressed as 

V = 2V ihl) + V {hh) + V {ir \ (32) 
whose evaluation is presented in the rest of this section. 



4.1 The heavy-light subsystem. 

As for the problem of evaluation the bare proper vertex V^ hl ' of composite operator (qQ) with a 
massless quark field q and the effective static heavy quark field Q, we can easily see that the result 
does not depend on the Dirac structure of the vertex. For this reason, in our calculations we have 
used the result of 0|, where this vertex was calculated to two-loop order in the Feynman gauge 



(a = 1) with the use of algorithm developed in [15 

CbQ>, 



T/(«) 
v 2,2 



*1,1 
1 



*1,0 



0. 



(33) 



Cb(-Cb — Ca) 



v (hi) 

V 2A 



-C B (C fl (l-4C(2))-C7 x (l-C(2))). 



Then from the relations ( PTj ) one can calculate the coefficients Z n ^, which determine the two-loop 
anomalous dimension for the subset of the heavy-light graphs 



7( W ) 



Cl(4-16C(2))-C B C A (4-4C(2)). 



(34) 



It is worth to note that this expression was calculated for antisymmetric baryonic color configuration 
q t Q j Q k eijk unlike the case of colour-singlet c^Q^S* mesonic configuration. The expression for the 
latter case can be obtained by substitution of Cb — > Cp, which reconstructs the required results, 



as it was checkeds by authors of [1 I 



4.2 The heavy-heavy subsystem. 

To evaluate this contribution we have used the results of [11], where the expression for the anomalous 
dimension of NRQCD mesonic vector current was presented 

■ j d In Z j , . 

if = 27« + 7^ = ^f (35) 



-C F (2C F + 3C„)y (^) 2 + 0(a 3 ,). 



In our case we have the similar problem, but the different color structure. Thus, following [I4| 
we can consider the matching of QCD vector current with the antisymmetric color structure on the 
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NRQCD one. Unlike the meson case with the singlet-color structure, the QCD vector current with 
the antisymmetric color structure need not be conserved, so we allow for its renormalization. In 
terms of the on-shell matrix elements, the matching equation can be written down as[] 



J 2,QCDAj,QCD l QCD 



CqZ- 



2,NRQCD^j^ NRQ CD l NRQCD 



+ 0(v 2 ), 



(36) 



where Zj^qcd has the following expression JTT 



'J, QCD 



17C A + 3C B + 4(1 + N F )T F ) 
- A C F (-17C A + 3C F + 4(1 + N F )T F )) (j^)' 



+ 



(37) 



The anomalous dimension of NRQCD current, obtained in this way, may be used in the calcu- 
lations of anomalous dimensions for the baryonic currents with two heavy quarks, as it does not 
depend on the Dirac structure of the vertex. The contributions of different two-loop diagrams with 
the antisymmetric color structure of the vertex q l Q^Q k tijk in the notations of |T3| are shown in 

(2} 

Appendix. To obtain the anomalous dimension 7}^) of composite operator under consideration one 
has to perform the following steps: 

1) sum all of these contributions, including the one- loop term, multiplyed by the two- loop 
QCD on-shell wave function renormalization constant |22]| , ZJq CD and one-loop NRQCD-current 
renormalization constant, 

2) perform the one-loop renormalization of coupling and mass. 

1 (2\ (2} 

After these manipulations the coefficient at - multiplied by —4 will give us the sum 7^ +2jq . 

(2) 

For the two-loop anomalous dimension 7^ in the heavy-heavy subsystem, we find the following 
result 



!(hh) 



--C B ((-19 + Qn 2 )C A 



4{tt 2 C b + 2N f T f )). 



(38) 



4.3 The light-heavy-heavy irreducible vertex 



In this case one needs to calculate the three-quark irreducible vertex Vq T \ There are 8 diagrams 
in the two- loop order. We have shown four of them in Fig. 1, the other four can be obtained by 
exchanging two heavy quark legs. We set the heavy quarks off shell in order to avoid any infrared 
singularities. Using the partial fractioning of the integrand in momentum integrals and recurrence- 
we arrive to the following expressions for the diagrams depicted on Fig.l 

7T 2 1 216 + 35tt 2 - 48^ (2) (1) - 96^ (2) (2) 
~3e + 36 

(40) 

r 2, 



relations of [15 



V 



(ir)[l] 



V r 



(ir)[2] 



2-Cl 



0. 



Air 



£-1 



(39) 



(ir)[3] 



V 



(tr)[4] 







2- CI 



2. a 



Air 
Air 



1 

~ e 

1 2 
+ - 



-4- 



4 



3vr 2 



(41) 
(42) 



4 Since the matching coefficient contains only short-distance effects, the matching can be done by comparing the 
matrix elements of these currents over a free quark-antiquark pair of the on-shell quarks at a small relative velocity. 
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where ^ n ^(z) = d n i\){z) / 'dz n , ip(z) = T'(z)/T(z) and the factor of 2 accounts for the contributions 
of the remaining four reflected diagrams not included in Fig. I. For the Z-factors and anomalous 
dimension, we obtain 

Z% ] = -3 • Cf , (43) 



2 



7$ = -4^=8-^1(1 + ^). (44) 
4.4 Anomalous dimension combined 

Now we are ready to calculate the anomalous dimension of baryonic currents with two heavy quarks. 
As we have already said above it does not depend on the Dirac structure of the current under 
consideration. Collecting the results for the heavy- light, heavy- heavy and irreducible light-heavy- 
heavy vertices, we find 

l { v ] = -^C B ((-13 + 30C(2))C A + 6(-2 + 6C(2))C B + 87V F T F ). (45) 

And, finally, to obtain the full two-loop result for the anomalous dimension one has to add the 
anomalous dimensions of heavy and light quarks. The result is 

if = ^(-48(-2 + 6C(2))C| + C A ((104- 240C(2))C B -101C F ) (46) 
-Q4C B N F T F + C F (-9C F + 52N F T F )). 

With this formula we are finishing our analytical calculations. 
Landing to the SU(3) group of QCD, we get 

7 « = -4, (47) 

m 254 152tt 2 20 , . 

7 (2) = + — N F w -194.909 + 2.222A^ F , (48) 

9 9 9 

which indicate a rather strong sensitivity of those currents to the choice of reference scale \i. 



5 Conclusion 

We have calculated the two-loop anomalous dimensions of NRQCD baryonic currents with two heavy 
quarks in the leading order in both the relative velocity of heavy quarks and the inverse heavy quark 
mass. It is shown, that the results do not depend on the Dirac structure of the currents and on 
the 75 prescription used in the calculations. These results will be useful for derivation of QCD sum 
rules for baryons with two heavy quarks in the same static approximation in both the leading and 
next-to-leading orders. We suppose to address this problem in the nearby future. 

This work is in part supported by the Russian Foundation for Basic Research, grants 96-02-18216 
and 96-15-96575. 
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6 Appendix 



In this appendix we present the generalization of expressions for the hard contributions to the 
diagrams of Fig. 1 of with the antisymmetric color structure of vertex evaluated at the threshold 



q 2 = Am 2 . Below you can find the coefficients of (a s /ii) 2 (e lE mQ/ '(Att/j, 2 ))' 
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Figure 1: The two- loop contribution to the light-heavy- heavy irreducible vertex with the spinor 
lines directed outside. 
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